Equilibrium statistical mechanics for incomplete nonextensive statistics 
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The incomplete nonextensive statistics in the canonical and microcanonical ensembles is explored 
in the general case and in a particular case for the ideal gas. By exact analytical results for the ideal 
gas it is shown that taking the thermodynamic limit, with z = q/{l~ q) being an extensive variable 
of state, the incomplete nonextensive statistics satisfies the requirements of equilibrium thermody- 
namics. The thermodynamical potential of the statistical ensemble is a homogeneous function of the 
first degree of the extensive variables of state. In this case, the incomplete nonextensive statistics is 
equivalent to the usual Tsallis statistics. If z is an intensive variable of state, i.e. the entropic index 
g is a universal constant, the requirements of the equilibrium thermodynamics are violated. 
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I. INTRODUCTION 

In modern physics there exist alternative theo- 
ries for the equilibrium statistical mechanics: the 
Boltzmann-Gibbs statistics [l|, 01, the Tsallis statis- 
tics [1, the Renyi statistics [1, Q based on Renyi 
entrop y M the incomplete nonextensive statis- 
tics 0, • All of them are obtained from the usual 
Boltzmann-Gibbs statistics due to the modification of 
the statistical entropy formula, the norm equation for 
a distribution function, and/or the way of taking en- 
semble averages. These statistics are defined at the 
equilibrium by the Jaynes principle of maximum en- 
tropy 111]. Therefore, they have to satisfy all require- 
ments of equilibrium thermodynamics and, moreover, 
must be in agreement with probability theory and gen- 
eral principles of physics. The incomplete nonexten- 
sive statistics is defined by the Tsallis entropy for- 
mula |3|] , the modified norm equation for a phase dis- 
tribution function and a modified expectation value 
for dynamical variables [1, [l3|- It differs from the 
Tsallis statistics only in the definition of the norm 
equation and ensemble averages. Nevertheless, there 
exist some other variants of incomplete nonextensive 
statistics which are defined either by an entropy dif- 
ferent from the Tsallis one [T2| - [l5j or by incomplete 
expectation values, but with the usual norm equation 
for a distribution function [l6| . 

The main aim of this Letter is to investigate the re- 
lation of the original incomplete nonextensive statis- 
tics d, [l3| to equilibrium thermodynamics, and to 
compare this with the one valid for the Tsallis statis- 
tics. A necessary and sufficient condition for fulfilling 
the zeroth law of thermodynamics, the principle of 
additivity, the Euler theorem, and the Gibbs-Duhem 
relation is that the thermodynamic potential of the 
statistical ensemble be a homogeneous function of the 
first degree of the extensive variables of state [13, [ll] • 
This is achieved in the thermodynamic limit. For 
finite systems the thermodynamical potential is an 
inhomogeneous function. The correct definition of 
the thermodynamic limit for the Tsallis nonextensive 
statistics was first considered in Botet et al. [19, uOl 



and in [l7|, [ij, l2l|. The problem of the zeroth law 
of thermodynamics for the incomplete nonextensive 
statistics was discussed in Refs. |2l| - l3G |. 

The main relations for the incomplete nonextensive 
statistics in finite systems for the canonical and mi- 
crocanonical ensembles were derived in our previous 
paper [61 . In the present Letter we discuss whether the 
incomplete nonextensive statistics can be thermody- 
namically consistent, i.e., whether the thermodynam- 
ical potential of the statistical ensemble is a homoge- 
neous function of the first degree of its extensive vari- 
ables of state. We concentrate on the thermodynamic 
limit, with z being an extensive variable of state. 

This Letter is organized as follows. In Section 2 we 
briefiy outline the microcanonical ensemble for the in- 
complete nonextensive statistics in general and in par- 
ticular for the ideal gas, and check the zeroth law of 
thermodynamics. The same procedure for the canon- 
ical ensemble is given in Section 3. 



II. MICROCANONICAL ENSEMBLE 

A. General formalism {E, V, z, N) 

The equilibrium statistical mechanics is an incom- 
plete nonextensive statistics if the equilibrium phase 
space distribution function satisfies conditions im- 
posed by the Tsallis' statistical entropy [3|], 



S = -k dT 



.9-9" 



= k{z + l) j dVg^il- Q^), (1) 

Here g satisfies an incomplete norm equation, and 
modified expectation values are taken as A [l3| , 



{A) 



dr£>« = j dTg—, (2) 
dVg'^A = / dVg^A. (3) 



Here j — g'^ — g"^^ is the phase distribution func- 
tion, fc is the Boltzmann constant, g G R is a real 



2 



parameter, q € [0, oo], and z is the variable of state, 
z = (?/(!-<?). 

Consider the microcanonical ensemble [E, V, z, N). 
Using the thermodynamical method based on the fun- 
damental equation of thermodynamics |l7| instead 
of the Jaynes principle [Hj, we obtain the equilib- 
rium phase distribution function and the statistical 
ht [6^ 

f = = W-'A{H -E)^ ga, (4) 
W ^ j dT = J A{H - E)dT, (5) 

D 

where H is the Hamiltonian function, gQ is the dis- 
tribution function of the Gibbs statistics, z is the 
thermodynamical variable of state [17[ and A(e) is 



the function distinct from zero only in the interval 
< e < AE, where it is equal to unity, D be- 
ing the region of phase space restricted by inequality 
E < H < E + AE. These two methods provide us 
with the same results for the distribution function(see 
Ref. [TtI and Refs. [H iH, HI). In the following the 
subscript G indicates the Gibbs statistics. The ther- 
modynamical potential of the microcanonical ensem- 
ble {E, V, z, N), the entropy S, has the form Q 



S = k{z + 1) 



1 



Sc^klnW, (6) 



where Sg is the usual Boltzmann-Gibbs entropy in the 
microcanonical ensemble. Since the distribution func- 
tion ^ and the entropy ^ have equilibrium values, 
we obtain the fundamental equation of therniodynam- 



TdS = dE+ pdV + Xdz - fidN, (H) = E, (7) 



where 



1 
T 

X 



dS_ 
'dE 



V,z,N 



z + 1 
zTg 



e , 



J, .dS\ TS TgSg 



= T 



^ -T 



^ E.V,N ^ + 1 



dS 



dV 



PG, 



E,z,N 



dS 
dN 



(8) 
(9) 
(10) 

(11) 



E,V,z 



where the definitions for the Gibbs measures can be 
found in Ref. 33i] . The fundamental equation of ther- 
modynamics ([7]) provides the first and second princi- 
ples of equilibrium thermodynamics 

5Q = TdS, 6Q = d{H) + pdV + Xdz - fj,dN, (12) 

where 6Q > is a heat transfer by the system to 
the environment during a quasistatic transition of the 
system from one equilibrium state to a nearby one. 

The zeroth law of thermodynamics, the Euler theo- 
rem and the Gibbs-Duhem relation for the incomplete 
nonextensive statistics in the microcanonical ensemble 



are valid whenever the thermodynamic potential S is 
a homogeneous function of the first degree of all ex- 
tensive variables of state For the homogeneous 
entropy S all functions of variables of state are either 
homogeneous functions of the first degree (extensive 
or additive of the first degree A = Ai + A2) 



A{E, XV, Xz, XN) 
A{E, XV, z, XN) 



XA{E,V, z,N) (z-extensive), (13) 
XA{E,V, z,N) {z-intensive) (14) 



or homogeneous functions of the zero degree (intensive 
or additive of the zero degree (j) — (pi ~ (f>2) 

(j){E,XV,Xz,XN) = (/)(£;, y,z,iV)(z-extensive), (15) 
(l){E,XV,z,XN) = </)(£;, F,z,iV)(z-intensive). (16) 

Here A is an appropriate scaling parameter; it can be 
A = \/E, X = l/N, X = \/V { or X = l/z for z 
extensive). Note that the extensive function A is ad- 
ditive, A = A1+A2, but the function is intensive, 
XA — XiAi = X2A2- The variables of state should 
satisfy the following conditions in equilibrium: 



E 
z 
z 



Ei+E2,N ^ Ni +N2,V 
z\ + Z2 (z-extensive), 
z\ = Z2 (z-intensive) . 



V1 + V2, 



(17) 
(18) 
(19) 



For extensive variables E,V,z,N the specific ratios 
XE, XV , XN and Xz are intensives. If the entropy 
5* is a homogeneous function of the first degree of all 
extensive variables of state , dH]) , then for the in- 
tensive functions in the microcanonical ensemble ([8]), 
(fTO]) . (fTT|) . we have 



T = Ti = T, 



p = Pi=P2, H = pi=^l2- (20) 



Therefore, the zeroth law of thermodynamics, T = 
Ti = T2, is satisfied whenever the thermodynamic po- 
tential S of the microcanonical ensemble is a homo- 
geneous function of the first degree of the extensive 
variables of state [13] ■ In this case, the entropy, its dif- 
ferential and the heat transfer are additive functions 
S = Si+ S2, dS = dSi + dS2 and 5Q = 5Qi + SQ2, 
respectively. The last equation follows from eq. ([H]). 
The homogeneous S satisfies the Euler theorem and 
the Gibbs-Duhem relation 

TS = (H) +pV + Xz - fxN (z-extensive), (21) 
SdT = Vdp + zdX - Ndfj, (22) 



and 



TS = {H)+pV-fiN (z-intensive), 
SdT = Vdp- Xdz - Ndfj,. 



(23) 
(24) 



For z being extensive, the Euler theorem (|21l) contains 
the term Xz and it is consistent with the fundamen- 
tal equation of thermodynamics ([7]) as well as with 
the Gibbs-Duhem relation (22). However, for z being 
intensive the Gibbs-Duhem relation (l24t is violated, 
and the Euler theorem (|23|) is not consistent with the 
fundamental equation of thermodynamics ([7]). Thus, 
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in this case the variable of state z can only be exten- 
sive. 

The Tsallis entropy ([T]), (jG]) may still become a 
homogeneous function in the thermodynamic limit. 
For z being extensive in this limit, functions of the 
variables of state are expanded in power series un- 
der the conditions N — oo, E oo, V oo, 
z — oo and e = E/N — const, v — V/N = const, 
z ^ z/N ^ const (a > 0). One determines 

A{E,V,z,N) = iV[a(e,w,z) +0(iV^")], (25) 
cj)iE,V,z,N) = (b{e,v,S) + 0{N-"), (26) 

where e = E/N is the specific energy, v — V/N is the 
specific volume, z — z/N is the specific z, and a = 
A/N is the specific A. For z being intensive, functions 
of the variables of state are expanded in power series 
under similar conditions as before, but with z = const. 

A{E,V,z,N) = N[aie,v,z) + 0{N-°')], (27) 
ct>{E,V,z,N) = <j>{e,v,z) + 0{N-"'). (28) 

One concludes that the entropy is a homogeneous 
function of the first degree in the thermodynamic 
limit. 



B. Ideal gas in microcanonical ensemble 

Let us investigate the thermodynamical proper- 
ties of the microcanonical ideal Maxwell-Boltzmann 
gas [6] in terms of the variable z — q/{l — q). The 
phase distribution function and the statistical 
weight (O for the ideal gas in the microcanonical en- 
semble can be found in Ref . ^33i] . The thermodynamic 
potential ^ is found to be 



S = k{z + l) (l - W- 



l/z 



(29) 



where Sg — k\nW is the Gibbs entropy. The tem- 
perature (O, the pressure (jlOp. the variable X ^ and 
the chemical potential pT|) are now 



T = To 

X =. kTa 



z + 1 ' 

z{W^/^ - 1) 



z + 1 



(30) 
flG, (31) 



where the definitions for the Gibbs measures can be 
found in Ref. [3l|. The heat capacity is given by 



Ce = 



k{z + l) 



kzTc 
E 



(32) 



Its general microcanonical definition can be found in 
Ref. [33i. In the Gibbs limit the entropy the 
temperature and the variable X can be written as 

S\z^±oo — So, T\z^±oo — Tq and X\z^±oc — 0, 
respectively. 

Let us investigate the thermodynamics of the ideal 
gas in a finite system. We divide the system into 
two parts, under the conditions P7)) -([T5 1) . Then from 



Eqs. (129j)-(130l) it follows that the Tsallis entropy for a 
finite system is nonadditive (nonextensive) function, 
S* 7^ 5*1 + 5*2 , and its proper temperature is noninten- 
sive function, T ^ Ti ^ Ti. Accordingly, the micro- 
canonical phase distribution function for this statistics 
does not factorizc, Q1Q2, where W ^ WiW2- Us- 
ing Eqs. ([29l) -([3T |) . we obtain the Euler theorem for 
the incomplete nonextensive statistics in the case of z 
being extensive 



TS + 0\^ = E+pV + Xz~fiN. 



O 



W 



Or 



kT, 



G 



z + l 



and in the case of z being intensive 



E + pV - fiN, 



(33) 
(34) 



(35) 



O 



= Og- kTaziW^^'' ~l--\nW], (36) 



where 
Og = E + kTcln 
+ kTGN 



ENIT { -N 



l-,p{N + l)-^^(^N 



(37) 



We conclude that for a finite ideal gas in the mi- 
crocanonical ensemble, using the incomplete nonex- 
tensive statistics, the zeroth law of thermodynamics, 
the principle of statistical independence and the Euler 
theorem are violated. They are polluted by finite-size 
effects. 



C. Microcanonical ideal gas in the 
thermodynamic limit with z being intensive 

Let us now consider the ideal gas in the thermody- 
namic limit with z being intensive. The phase distri- 
bution function g and the statistical weight W in the 
thermodynamic limit can be found in Ref. [33j . It is 
readily seen that the entropy ([29]) . the temperature 
and pressure ([50]) . the variable X and the chemical 
potential ((3T|) can be written as 



S = k{z + 1) (l-w-^/^) 



(38) 



and 
T 
X 



Tg 

kTG 



P = PG, 



Z + 1 ' 

z(w^/^ - 1) 



Nlnw 



(39) 

, M = /iG, (40) 



where the Gibbs measures in the thermodynamic limit 
can be found in Ref. [33]. Finally, the heat capacity 

(I32D is found to be Ce = k{z + l)u;-^/^ (l + |f )~^ 
In the Gibbs limit the incomplete nonextensive statis- 
tics quantities (p8)) - p0)) recover their Gibbs val- 
ues, the entropy S\z~^±oo = Sg, the temperature 
T\z^±oo = Tg, the variable X\z^±oo — and the 
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heat capacity Ce\z^±oo = Ce,g = 3iV"fc/2, respec- 
tively. Therefore, for z being intensive in the mi- 
crocanonical ensemble the Gibbs limit is commutative 
with the thermodynamic limit. 

Considering Eqs. pT| . (IT9)) . e = £i = £2 and 
V = vi = V2, in the thermodynamic hmit for z being 
intensive, we obtain that the entropy (j38p is nonad- 
ditive. The temperature p9p is nonintensive and the 
phase distribution function factorizes, g = QiQ2- The 
Euler theorem (j35p for the incomplete nonextensive 
statistics is violated 



= -kTGz[w 



N/z 



1 



N 



■ In-i 



(41) 



We conclude that the zeroth law of thermodynamics 
is not satisfied for the ideal gas in the thermodynamic 
limit with z being intensive. 



D. Microcanonical ideal gas in the 
thermodynamic limit with z being extensive 

Finally, we consider the thermodynamic limit with 
z being extensive. The phase distribution function g 
and the statistical weight W in the thermodynamic 
limit are calculated as in Ref. [s^]. The entropy (P^ . 
the temperature and pressure (1301) , the variable X and 
the chemical potential (j3ip can be written as 



5* = kzN 



(1 - 



and 



T = IG 

X = kTG 



P ^ 
1/5 _ 1 _ 



PG, 



fl = ^IG- 



(42) 

(43) 

(44) 



The Gibbs quantities, labelled by the index G, are 
given in Ref. [33|. The heat capacity ([32]) takes 

the form Ge = kiNw'^''^ (l + f )~^ In the Gibbs 
limit the incomplete nonextensive statistics quantities 
recover their Gibbs values, the entropy S\i^±ao = 
Sg, the temperature T|i_j.±oo — Tg, the variable 
±00 — and the heat capacity G e\z^±oo — 
Ge,g = 3iVfc/2, respectively. Therefore, for z being 
extensive the Gibbs limit commutes with the thermo- 
dynamic limit. Let us comment that the results P^- 
for the incomplete nonextensive statistics coin- 
cide with res pec tive results for the Tsallis statistics 
given in Ref. [l7| . 

Considering Eqs. ([T7|) . (fTS)) and e — ei — 82, 
V — vi = V2, and z = zi = Z2, in the thermody- 
namic limit for z being extensive, we obtain that the 
entropy (j42|) is additive, the variable X (j44|) and the 
temperature (|43p are intensive and the phase distri- 
bution function factorizes, g — gig2- Moreover, the 
Euler theorem ([55)) is valid, = 0. Thus, for the 
incomplete nonextensive statistics in the thermody- 
namic limit with z being extensive the entropy is a 
homogeneous function of the first degree and the ze- 
roth law of thermodynamics, the Euler theorem and 
the principle of statistical independence are valid. 



III. CANONICAL ENSEMBLE 

A. General formalism (T,V,z,N) 

Let us investigate the canonical ensemble 
(T, y,z,iV). To obtain the equilibrium phase 
distribution function, we use the thermodynamical 
method explored in [13, [ll] instead of the Jaynes 
principle 



1 



(z + l)2 kT 



(45) 



where z = q/{l — q) and A = A(T, V, z, N) is the norm 
function which is the solution of Eq. For details, 
see Ref. @. These two methods provide us with the 
same results for the distribution function (see Ref. [l^ 
and Refs. [HI, iH, [H). The expectation value (A) of 
a dynamical variable A is determined by Eq. (|31) with 
the distribution function (j45p . The entropy and the 
free energy are expressed through the norm functions 
A and (H) as 



S 
F 



^((g)-A) 
T{z + 1) ' 

{H) -TS^ 



(H) + zA 
z + 1 



(46) 

(47) 



Equation ()47p is the Legendre transform of energy 
with respect to the entropy of the system. The funda- 
mental equation of thermodynamics in the canonical 
ensemble is the same as Eq. ([7]) but with 6\ 



T,z,N 



dN 



T,V,z 



X ^ kT dVg^ 



-2 + 1 ^ , -J— 

1 - p-'+i H gi^+i Inp^+i 



(48) 
(49) 

(50) 



where the following properties of the Hamiltonian 
function, {dH/dT)v,z,N — {dH/dz)Ty,N — 0, were 
assumed. It provides the first and the second laws 
of thermodynamics (|12p . The thermodynamical rela- 
tions for S, p, X and ^ in the canonical ensemble can 
be found in Ref. Ts]. Note that the thermodynami- 
cal relations for the c^uantum incomplete nonextensive 
statistics are the same as for the classical one given 
above. 

The zeroth law of thermodynamics, the Euler theo- 
rem and the Gibbs-Duhem relation for the incomplete 
nonextensive statistics in the canonical ensemble are 
valid if the thermodynamic potential is a homoge- 
neous function of the first degree of all extensive vari- 
ables of state. The general proof of this statement is 
the same as the one given in Refs. [l^, H^. Moreover, 
the definition of the thermodynamical limit and the 
proof of the extensivity for the variable of state z in 
the canonical ensemble are the same as for the Renyi 
statistics 1331. 
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B. Ideal gas in the canonical ensemble 

Let us rewrite the main relations for the ideal gas 
in the canonical ensemble @ in terms of the variable 
of state z = q/{l — q). The phase space distribution 
function (|45|) and the norm function A can be written 
as 



N pI_ 
2m 



[z + lf 



kT 



A 
and 



kT 



(z + 1)2(5-1) 



(51) 



(52) 



0'z,N 



0.z,N 



r(z + i) 



r(-z- |7V) 



z > 0, (53) 



((iTfF)*^r(-z) 



z < 0, 



(54) 



where g is the spin-isospin degeneracy factor and z < 
-|iV in Eq. JSl]). The partition function Zq for the 
ideal gas was defined in Ref. fl§\ . The thermodynamic 
potential can be written as 



F 



kT 



1 



3N 

2(2+1) 



zA 



3N 



2z {z + l)kT 



(1 



(55) 



The energy, the entropy (j46|) and the pressure (148 
can be written as 



(H) 



2 z i + i-4^ 



2 z + 1 



-N - 

2^* 



z A 
z + 1 kT 



P 



^ 2 z+1 

V z 



^ ^ 2 z+1 



(56) 
(57) 



The chemical potential is found to be 

(H) 



N{z + 1 + ^N) 



N 

where A,, 



ln(l 



A 



{z + lfkT' 
d{\naz^N)/dN and 



kT' 



,(59) 



3 

^ - -7/;(z + 1 + -A^) - ln( 



3 
2A„ 



(^^)'^>0'(60) 
),z<0. (61) 



Here ip{z) is the psi- function, z < — |Af in Eq. (|5l]) . 
The chemical potential /^g for the Gibbs statistics was 
given in Ref. |33j . The function X can be written as 



(H) ^ + Ii£tv zA 

^(^ + 1) 1 + (z + iy 



x„ 



2{H) 
3 Af 
37V(z- 1) 
2z(z + 1) 



+ In 1 



A 



(z + l)2 kT^ 



(62) 



3A^ 

V'(z + 1 + — )+i^(z + l), z>0, (63) 

3Ar(z- 1) 
2z(z + 1) ~ 

37V 3 
- V-l-^- — )+i^(-^), ^<-2^' 

where ATa = 9(lna^jv)/9z. The heat capacity is now 
CvzN = {H){l + ^)~^/T. Its definition in the canon- 
ical ensemble can be found in Ref. [33*1 . In the Gibbs 
limit z — >■ ±oo,A^ = const, V = const the incomplete 
nonextensive statistics recovers the usual Boltzmann- 
Gibbs statistics. We conclude that for the finite ideal 
gas in the canonical ensemble the thermodynamical 
potential (j55l) is an inhomogeneous function. There- 
fore, the zeroth law of thermodynamics, the Euler the- 
orem and the principle of additivity are not satisfied. 



C. Canonical ideal gas in the thermodynamic 
limit with z being extensive 

Let us consider the thermodynamical limit with z 
being extensive, A^ — oo, V — > oo, z — ±oo and 
V = V/N = const, z = z/N = const. The phase 
distribution function (1511) and the norm function A 
((52|) take the form 



/ 
A 



= Q = 



1^1=1 2m 



IN 



kT 



IN 



(65) 



kTlN{B -\),B = {Zgc""^) (1 + ^),(66) 



N 



3 

^+2i 



-(z+f)Ar 



(67) 



The partition function Zq, Zq can be found in 
Ref. [33. The thermodynamic potential (1551) in the 
thermodynamic limit with z being extensive takes the 
form 



F = kTzN{B - 1) A. 



(68) 



Therefore, in the thermodynamic limit the norm func- 
tion A coincides with the free energy, because in the 
series (j68p only the term proportional to the first 
power of N is kept. The thermodynamical potential 
((55)) is a homogeneous function of first degree since z 
is a intensive variable. The energy (j56p . the entropy 
(j57p and the pressure (j58p in the thermodynamic limit 
can be written as 



S 



kNz 1 



B 



1 + ^ 

' 2z 



-NkT- 



B 



1 + 51 



kT B 



(69) 
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The energy ([70)) and the entropy ([69)) are extensive 
functions but the pressure p is intensive. The variable 
of state X ([62)) and the chemical potential ([59)) can 
be written as 



X 



kT 



B 



1 + A 



In. 



B 



1 + 



-kT- 



B 



1 + 



2 1 + 



B 



_3_ 

2z 



fcT 



(71) 
(72) 



The chemical potential ^.q in the thermodynamic 
limit can be found in Ref. j33:]. The heat capacity 
is now CvzN = {H){1 + ^)"VT. In the Gibbs limit 
z — > ±oo and N — const the incomplete nonextensive 
statistics recovers the usual Boltzmann-Gibbs statis- 
tics. For example, the measures ([68)) . ([69)) take the 
form F\i^±oo = -Fb, S\i^±oo = So, respectively. 
This way in the incomplete nonextensive statistics, 
with z being extensive, the Gibbs limit commutes with 
the thermodynamic limit. This conclusion differs from 
the results found in Ref. [3^ . 

Dividing the system into two parts under the con- 
ditions T = Ti = T2, V = vi = V2 and z — zi = Z2, 
we obtain that the Tsallis entropy ([M)) is a homo- 
geneous function of first degree (extensive) and it is 
additive, S = Si + 82- Therefore, for the ideal gas 
in the canonical ensemble in the thermodynamic limit 
with z extensive the zeroth law of thermodynamics, 
T = Ti = T2, and the Euler theorem ^ are valid. 
Let us note that in the thermodynamical limit with z 
extensive all relations for the ideal gas in the canonical 
ensemble (l55t - ([7^ completely coincide with the ones 
of the original Tsallis statistics given in [l^ . Thus the 
incomplete nonextensive statistics satisfies all require- 
ments of equilibrium thermodynamics in the canonical 
ensemble in the thermodynamical limit with z exten- 
sive. The functions of the variables of state are ho- 
mogeneous functions of first degree, for extensive, or 
homogeneous functions of the zero degree, for intensive 
variables. In particular, the temperature is an inten- 
sive variable and thus provides implementation of the 
zeroth law of thermodynamics. 



D. Canonical ideal gas in thermodynamical 
limit with z being intensive 



We inspect an ideal gas in the thermodynamical 
limit with z being intensive: N ^ 00, V ^ 00 and v = 
V/N = const, z = const. For the sake of convenience 
let us consider only the case z > 0. The phase space 
distribution function ([5T)) and the functions ([5^ -([5 ^ 
take the form 



{z + iy 



kT 



A = -NkTZ. 



-2/3-1 



(73) 
(74) 



where 



az,N 



r(g -H)e^+i+i^/72^ 



(75) 



N 



The thermodynamic potential ()55p in the thermody- 
namic limit with z being intensive takes the form 
F ^ A. Therefore, the norm function A coincides 
with the free energy. In the series expansion of F 
only the leading term proportional to the first power 
of N is kept. The energy ([55)) . the entropy ([57)) and 
the pressure ([58)) can be written as 



{H) 
0{N°) 

P 



A, 



k{z + l)[l- 



z + 1 



-Z. 



G 



V ^ 3 V 



(76) 
(77) 

(78) 



Therefore, the entropy S — 0. The function of state 
X ([62)) and the chemical potential ([59)) can be writ- 
ten as AT = and /Lt = | ^ . The heat capacity is now 

CvzN = O(iVO), where 0{N°) = kZ~'^^^e-\ In the 
Gibbs limit z — ^ ±00 and N = const the distribution 
function ()73l) and the norm function ([7^ do not re- 
semble the distribution function and the free energy 
of the Gibbs statistics. The variables F, {H) , S, p, X, fi 
and CvzN also do not resemble their Gibbs values. 
Thus, for the incomplete nonextensive statistics with 
z being intensive the Gibbs limit does not commute 
with the thermodynamic limit as in Ref. [33|. This 
means that the thermodynamic limit with z being in- 
tensive is erroneous for both the incomplete nonexten- 
sive statistics and the Tsallis one [33 ]. 

The Euler theorem ([23)) is valid if for extensive func- 
tions we keep the terms proportional to the first power 
of N and for intensive ones we keep the terms pro- 
portional to N'^. In the thermodynamic limit with z 
being intensive, the entropy S is intensive and it is a 
homogeneous function of the zeroth degree. This con- 
tradicts the requirements of the equilibrium thermo- 
dynamics. Moreover, the thermodynamic limit with z 
being intensive is not consistent with the Gibbs limit. 
Therefore, the ideal gas results with z being intensive 
completely contradict the equilibrium thermodynam- 



IV. CONCLUSIONS 

In this Letter, we have examined the incomplete 
nonextensive statistics in the canonical and micro- 
canonical ensembles in the general case and in the 
particular case of the classical ideal gas. The exact 
analytical results for the ideal gas were obtained for 
both finite systems and in the thermodynamic limit. 
The ideal gas in the thermodynamic limit was stud- 
ied in two cases: the variable z being extensive or 
intensive. Connections between incomplete nonexten- 
sive statistics and equilibrium thermodynamics were 
analyzed. 
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To this end, we summarize our main results. The 
phase distribution functions for the canonical and mi- 
crocanonical ensembles were derived by the method 
based on the fundamental equation of thermodynam- 
ics. This method gives the same result as starting 
from the Jaynes principle. The connection of the 
incomplete nonextensive statistics with equilibrium 
thermodynamics was tested. The main thermody- 
namical relations were obtained from ensemble av- 
erages by using phase space distribution functions. 
It is shown that the zeroth law of thermodynamics, 
the principle of additivity, the Euler theorem and 
the Gibbs-Duhem relation are valid only if the cor- 
responding thermodynamical potential is a homoge- 
neous function of the first degree of all extensive vari- 
ables of state. By analytical results for the ideal gas it 
was revealed that the thermodynamical potential can 
be a homogenous function of the first degree only in 
the thermodynamic limit, whenever the entropic in- 
dex z is an extensive variable of state. Moreover, it 
was found that the variable z must be extensive in 
order the Euler theorem be consistent with the funda- 
mental equation of thermodynamics and the Gibbs- 



Duhem relation. 

On the whole, our analytical results for the ideal gas 
indicate that the incomplete nonextensive statistics in 
the canonical and microcanonical ensembles is ther- 
modynamically consistent whenever it coincides with 
the usual Tsallis statistics in the thermodynamic limit 
with z being an extensive variable of state. Therefore, 
we conclude that the incomplete nonextensive statis- 
tics in the equilibrium statistical mechanics duplicates 
the thermodynamic relations of the familiar Tsallis 
statistics and leads to no new result in this respect. 
This conclusion relies on the particular approach that 
the extra parameter q in the entropy formula is treated 
as related to a variable of state, z — q/{l — q). 
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